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Abstract. We calculate the power emitted in scalar modes for a binary systems, 
including binary pulsars, with a conformal coupling to the most general Galileon effec- 
tive field theory by considering perturbations around a static, spherical background. 
While this method is effective for calculating the power in the cubic Galileon case, here 
we find that if the quartic or quintic Galileon dominate, for realistic pulsar systems 
the classical perturbative expansion about spherically symmetric backgrounds breaks 
down (although the quantum effective theory is well-defined). The basic reason is 
that the equations of motion for the fluctuations are then effectively one dimensional. 
This leads to many multipoles radiating with equal strength, as opposed to the nor- 
mal Minkowski spacetime and cubic Galileon cases, where increasing multipoles are 
suppressed by increasing powers of the orbital velocity. We consider two cases where 
perturbation theory gives trust- worthy results: (1) when there is a large hierarchy be- 
tween the masses of two orbiting objects, and (2) when we choose scales such that the 
quartic Galileon only begins to dominate at distances smaller than the inverse pulsar 
frequency. Implications for future calculations with the full Galileon that account for 
the Vainshtein mechanism are considered. 
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1 Introduction 



Attempts to explain the observation of a small but nonzero vacuum energy in a dynam- 
ical way generically introduce new light scalar degrees of freedom, [1, 2]. Frequently 
these scalars are themselves radiatively unstable. For example, in massive gravity, the 
new helicity-0 polarization state of the graviton acts as a light scalar. As such light 
degrees of freedom have not been detected in precision tests of gravity, the phenomeno- 
logical viability of these theories rests on the existence of screening mechanisms that 
hide the scalars on small scales where gravity has been well tested. Therefore it is 
imperative to understand in detail how screening mechanisms work. In this work we 
will study how one such screening mechanism, the Vainshtein mechanism [3], operates 
around dynamical, non-relativistic sources. Vainshtein screening works by making the 
scalars strongly coupled to themselves near compact sources, essentially suppressing 
their coupling to matter. 

While studying the Dvali-Gabadadze-Porrati (DGP) model of 'soft' massive grav- 
ity [4], a class of scalar field models was discovered that exhibit the Vainshtein mecha- 
nism. These fields posses the galilean symmetry 7r — > ir + c + v /1 x fl and are thus referred 
to as 'Galileons' [5]. They arise naturally in the decoupling limit of theories of massive 
gravity, both 'soft' such as DGP [6], and 'hard' such as massive gravity from auxiliary 
extra dimensions, [7-9], ghost free dRGT massive gravity [10, 11], or New Massive 
Gravity in three dimensions, [12, 13]. They can also be considered as an effective field 
theory in their own right, which is the approach we will take in this work. Vainshtein 
screening for Galileons around static, spherically symmetric is well understood and is 
quite effective [5, 6]. There has recently been a lot of work studying the Vainshtein 
mechanism [14-34]. 

However, recently it has been established that within its simplest realization (i.e. 
in the cubic Galileon), the Vainshtein mechanism is slightly less effective around dy- 
namic sources, because the orbital period introduces a new large length scale flp 1 [35]. 
There it was shown that the Vainshtein suppression to the gravitational radiation is 
(fipr*) -3 / 2 , instead of the naive expectation from static sources (f/r*) 3 / 2 , where f is 
the size of the system and r* is the Vainshtein radius. Nevertheless, the Vainshtein 
mechanism in that case ends up being still powerful enough to evade any constraints 
from pulsar systems. A similar problem was also recently considered in [36]. 

In this work we extend the results to the power emitted by binary pulsars in the 
presence of all four non-trivial Galileon terms possible in four dimensions. We consider 
a conformal 7rT coupling to matter. We find that, surprisingly, the naive perturbation 
method used for instance in [37] fails to work in this case. This is despite the fact that 
the perturbative expansion is valid for the cubic Galileon case [35] as we show explicitly 
here in Appendix E. The difference is that when the quartic Galileon dominates, the 
effective one dimensional metric seen by the fluctuations means modes radiate with an 
effective frequency 



Multipoles radiate as long as their associated effective frequency uj£ is real. The crucial 
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point is that there is no distinction between modes of different £ when t <C wr^, so 
each multipole radiate with equal strength. Since Qpr+ ^> f (for the Hulse- Taylor 
pulsar £lpr* ~ 10 6 , when considering a strong coupling scale A ~ (lOOOkm)" 1 ), a 
huge number of modes radiate, each with a comparable power. Even more multipoles 
radiate when we consider the power emitted by higher harmonics with frequency nflp 
for integer n. 

We interpret this result as a failure of perturbation theory, and explicitly demon- 
strate that the perturbation series breaks down. While we are able to find some 
regimes where perturbation theory is valid and obtain sensible results, the mere fact 
that perturbation theory breaks down in this situation questions the intuition one can 
really infer from the static Vainshtein mechanism when dealing with more complicated 
systems. 

We stress that this breakdown of perturbation theory is not a quantum strong 
coupling problem. The quantum low energy effective theory for the Galileon remains 
under control, rather this is the breakdown of the description of the classical field con- 
figuration for the Galileon around the pulsar as a nonlinear spherical profile plus small 
time-dependent non-spherical perturbations. It is thus a failure of the approximation 
method used to calculate radiation and not of the effective field theory per se. 

These results suggest that there is additional Vainshtein suppression on top of 
the naive expectation from the static, spherically symmetric case. We obtain an es- 
timate for the power by introducing a cutoff I in the number of multipoles that we 
sum to determine the power. Physically this is motivated by the fact that the higher 
order multipoles are more strongly angular dependent and so more sensitive to the 
non-spherical nature of the system. We expect that multipoles at higher I will be- 
come nonlinear and see Vainshtein screening on top of the Vainshtein screening from 
the background. With this procedure we find that the total power obtained is more 
Vainshtein suppressed ~ r~ 3 than any one mode ~ r~ 2 . 

The rest of this paper as organized as follows. In Section 2 we derive the pertur- 
bation equations around a static, spherically symmetric background and review how 
to use the Feynman propagator to derive the power emitted. In Section 3 we apply 
the naive perturbative methods of [35, 37] to compute the power in a binary pulsar 
system, and find that in that case the resulting power is formally divergent. We then 
check whether perturbation theory is under control by explicitly constructing the first 
and second order solutions using the retarded Green function and comparing them. 
We find that, unlike the cubic Galileon case [35] the orbital velocity small parameter 
v is not enough to ensure the convergence of the perturbation series, and one needs an 
additional small parameter to have a controlled expansion around a spherically sym- 
metric background. This can take the form of a hierarchy between the two masses in 
the binary system, or between the strong coupling scale for the cubic Galileon and the 
strong coupling scale for the quartic Galileon. We consider these two cases in sections 
4 and 5 respectively as a check that the Galileon theory is well behaved in regimes 
where we can trust perturbation theory. We conclude by considering possible methods 
to calculate the power emission and the implications of this result for other studies of 
the Vainshtein mechanism. 
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2 Perturbations in a Galileon Theory 



We consider a Galileon theory around Minkowski spacetime. These theories can arise 
from fully covariant theories in a number of different ways [38, 39]. For example one 
can start with the action for ghost-free massive gravity [11] 

S = ^j dW^ (r - ^U(g, H)j , (2.1) 

where the tensor H^ u is defined in terms of the metric and the four Stiickelberg fields 

<p a as 

9^ = d^ a d„(j> b Va b + . (2.2) 

The ghost-free potential U(g, H) is fixed up to two free parameters (in addition to the 
mass parameter m). The Galileon effective field theory then arises in the decoupling 
limit Mpi — > oo, m — >■ keeping the scale A = (m 2 Mpi) 1 ^ 3 fixed, where the helicity-0 
mode plays the role of the Galileon. 

Regardless of the full covariant completion, we start with a Galileon scalar field 
Ti in Minkowski spacetime conformally coupled to matter [5] and consider all possi- 
ble interactions that respect the Galileon properties in four dimensions (ignoring the 
tadpole), 

where (Sh)^ = — ^\3h^ u + • ■ ■ is the Lichnerowicz operator, and T is the trace of the 
stress-energy tensor. As we can see, the helicity-2 and -0 modes decouple in this case 1 . 
Here the are the Galileon interactions in four dimensions 

C 2 = (<9tt) 2 (2.4) 
C 3 ee (<9tt) 2 [II] (2.5) 

£ 4 ^(avr) 2 ([n] 2 -[n 2 ]) (2.6) 
£ 5 ee (dixf ([n] 3 - 3[n] [n 2 ] + 2 [n 3 ]) (2.7) 



where IT^ = d^d^ and II" = n^II^ . . . n" n . Square brackets [A] denote the trace 
of the tensor A with respect to the Minkowski metric [A] = rj^A^. 
The equations of motion for n and h^ v are then 

(Sh)^ = -±-T, v , (2.8) 



3 v^U d 8Ci 1 



4 ^ A 3( i_2) Sir 2Mpi 



T. (2.9) 



1 We point out however that in Ghost-free Massive Gravity the helicity-2 and -0 modes do not 
fully decouple in the 'decoupling limit' unless C5 = 0. Furthermore that theory also includes other 
non-conformal couplings to matter which we ignore here. See Rcfs. [20, 29] for effects arising from 
these non-conformal coupling to matter. 
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The coefficients Cj are dimensionless, and so far arbitrary, although 02,03,04 must be 
positive for the stability of the theory (see Ref. [5]). Without loss of generality, we can 
absorb C2 into the definition of ir and C3 into that of A (i.e. we can set Ci = 1 and 
C3 = 1/3) 2 . The Aj are the scales associated with each of the Galileon interactions. 
These scales are typically assumed to be of the same order A and in theories of ghost- 
free massive gravity this scale is related to the mass of the graviton by A ~ (rr^Mpi) 1 / 3 . 
Current bounds on the graviton mass typically require m not to be too much larger 
than m ~ H ~ 10 -33 eV, in which case A ~ (1000 km) -1 . However, for a generic 
Galileon theory, one could potentially consider these scales as being different. The 
non-renormalization theorem present in Galileon theories allows for such a hierarchy 
without fine-tuning issues, [5, 40]. The notation used here is similar to that in Ref. [35], 
after setting A 3 = A, c 2 = 1, c 3 = 1/3. For simplicity we also set c 4 = 1 in what follows. 
We leave C5 arbitrary because it can be of either sign. 

In the rest of this paper, we will assume that a Vainshtein mechanism does oc- 
cur and that at short enough distances the interactions (2.5) or (2.6) dominate over 
the standard kinetic term (2.4). This depends on the relation between the different 
coefficients in the theory, which we assume to be a given fact. 

If A4 < A3 then £3 never has any effect because it will only become relevant 
at energies where it is already dominated by £4. If on the other hand A3 < A4, 
then the interactions £3 can dominate for a little while before being taken over by 
£4 at short enough distances. As we will see, it will also be convenient to take A 5 > 
(A4/ A.3) 1//3 A4. In the first part of this paper we will have in mind the situation A4 < 
A3 ~ (1000km) -1 . Notice that for a spherically symmetric configuration, the quintic 
interactions A5 vanish, and so these interactions are only relevant at the perturbed 
level (however as we shall see, even at that level, they simply correspond to a rescaling 
of some parameters) . 

Our basic philosophy for computing the power emission is to perform a back- 
ground+perturbation split in the Galileon where the background is static and spher- 
ically symmetric and the deviations from spherical symmetry is captured by the per- 
turbations. In effect this decomposition assumes that the majority of the Vainshtein 
screening comes from the monopole moment of the binary system. More precisely we 
split the field ir and the source as 



Physically this split is suggested from standard Effective Field Theory consider- 
ations, where we expect that the physics responsible for the radiation should arise at 
the (energy) scale Qp. Since this is scale is much smaller than the scale associated 
with the size of the system f -1 (typical distance between the two objects), where the 

2 Strictly speaking redefining c-i corresponds to changing the coupling to external matter. Of course 
by making that coupling small we reduce the amount of radiation into the Galileon, but in their natural 
realizations the coupling to matter is of order 1 and one is bound to rely on the Vainshtein mechanism 
to hide this scalar field. 
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spherical symmetry is broken, we expect that spherical background should be a good 
approximation when computing the power, barring some unusual circumstances. How- 
ever we emphasize that we are free to choose any background+perturbation split so 
long as the resulting perturbative expansion remains under control. 

The rest of this section is organized as follows. In the next two subsections we will 
solve for the background field it and derive the equations of motion for the fluctuations 
<fP-\ Then we will review how to compute the power using the Feynman propagator 
constructed from the fluctuations. 

2.1 Static and Spherically Symmetric Background 

Assuming a point source T M „ = -M tot 5 (3) (x)5g5°, where M tot = M 1 + M 2 is the total 
mass of the binary system, the background solution for 7r is spherically symmetric. 
Using the notation Vir(r) = rE(r), the background field equation for 7r takes the 
simple algebraic form 

E\ 2 fE\ 2 2 (E\ 3 _ 1 M tot 1 
7j + 3AfV7j + AfW ~12^Mp^' (2 ' 12) 

The quintic Galileon does not affect the background configuration, [5]. This is because 
the fc-th Galileon term is a topological invariant in dimensions smaller than k — 1, and 
since the system is static it is effectively three dimensional. One therefore has three 
branches of solution. We focus here on the 'normal' branch, which smoothly connects 
a free (weakly interacting) field E ~ 1/r 2 at spatial infinity r — > oo to a strongly 
interacting field at short distance scales, so as to achieve the Vainshtein mechanism. 

The source has two Vainshtein radii, r* )3 and r*^ associated to the two interaction 
scales A 3 and A 4 , 



( 



M tot \ 1/3 1 



r * s UsU M (213) 



*,4 



\16Mpi/ V A 4/ A 4 



(2.14) 



These two radii define three regimes in space, where £2,3,4 dominate in turn. More 
precisely, E(r), which is just the radial derivative of the background solution, is given 
by 

Mtot/Mpil 

E(r » r* j3 ) = 12?r ~ 2 (2-15) 

W( ^ ^ \ (Mtot/M P1 ) 1/2 A 3/2 , 01R , 

E(r* A < r < r^ 3 ) = == A 3 7 (2.16) 

Vo7rr 
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2.2 Equations of motion for perturbations 

The perturbed stress energy tensor 5T^ U which encodes the time dependent dynamics, 
which for slowly moving sources is 



ST.? = - 



Mi6 {3 \x - Xi(t)) - M tot 6 3 



x) 



•i=l,2 



y^ 



(2.18) 



where Mj is the mass of each companion 3 . The quadratic lagrangian for <ff& is 
£ = -\z^{x)d^d v ^ + -4^5T 



(2.19) 



pi 



l -Z tt {r){d t <P^ f - \z„{r) (d r ^) 2 - (Vn0 (1) ) 2 + -|^T (2.20) 



pi 



where 



Z rr (r) 



Zttir) 



4 E(r) 



3A 3 r 



J_d_ 

3r 2 dr 



r 3 I 1 



Y (\ 1 d 

2r dr 



r 1 + 



6 E(r) 2 

2 £(r) 
A| r 

4 £(r) 
3A 3 , r 



18 £(r) 2 
Af r 2 

6 E(r) 
f Af^ 2 " 



24£(r) 



(2.21) 
(2.22) 
(2.23) 



For notational consistency we include here A 5 . Recall however that A 5 does not con- 
tribute at the background level, and in particular there is no associated Vainshtein 
radius r* ; 5. Its only effect is to redress the time derivative pieces of the action. 

As we will see the quadratic action will only be sufficient to compute the power 
when there is a large hierarchy between A3 and A4 or between M\ and Mi. However 
for now we continue on without making any assumptions. The quadratic action gives 
rise to the equations of motion 



1 



5T. 



V6M Fl 

where □ is the modified d'Alembertian defined as 

□ = -d t (Z tt d t ) + d r (Z rr d r ) + A^Z nn V 



(2.24) 



(2.25) 



where is the Laplacian on a unit 2-sphere. It is useful to consider the form of □ in 
the different regions. 

• £2 region: r ^> 7^3 

□0« = D0M , (2.26) 



3 See Rcfs. [35, 36] for relativistic corrections to this expression. 
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i.e. far from the source the field is weakly coupled and perturbations are free as 
required. 

• £ 3 region: 7^4 <r< r* 3 
In this case, so long as A 5 ^> (A 4 /A 3 ) 1 / 3 A 4 the equations of motion reduce to the 
normal cubic Galileon equation of motion 



The case of small A5 is considered in Appendix B. 
• £4 region: r <C r* 4 



W) - 128 x 31/3 fM'f^ 2 
" 7T 2 / 3 U 3 y V r 



+ a r V + ^v 2 , 

C r r *,4 



(2.28) 



where the speed of sound of the radial fluctuations c r is given by 



4 A 



12 \ -1/2 



and the coefficient ka is given by 



6 1 



7r 2 / d / 27 /A 3 \ 
fe =i728^( 1 -TUj I' (2 ' 30) 

Note that in the £4 region the fluctuations effectively see a one dimensional metric 
ds 2 = —Z fJiU dx tJ- dx u oc — dt 2 + dr 2 + r 2 df2 2 , where crucially the angular part of the 
metric is multiplied by the constant r 2 instead of the normal factor of r 2 . 

Note that the second term is order 1 so long as A5 > (A4 / /A3) 1 / 3 A4. This is the 
same condition that we found above for the £5 contribution to be negligible in the £3 
region. The effect of £ 5 here is to decrease the sound speed, c r ~ (A5A3/A4 2 ). This 
case is considered in Appendix B. 

The stability of these theories was studied in Ref. [5], so as long as we take our 
coefficients to satisfy the conditions of Ref. [5] perturbations are guaranteed to be 
stable about the spherically symmetric configuration. 

2.3 Computing the power using the effective action 

Following [37] we compute the power in the binary pulsar system by looking at the 
imaginary part of the effective action 4 obtained by integrating out the fluctuations <p. 
We start with the quadratic action for the perturbations 

= f d 4 x^ (J-Z^d^^ -J—tfT) , (2.31) 



4 This method differs slightly from that followed in [36], but both strategies are valid and are 
ultimately equivalent. 
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where the are the trajectories of the two objects. We then integrate out the field 
leaving us with an effective action 

S cH = f d A x C M + jt^tjt [ d 4 xdV 'ST{x)G F {x,x')ST(x') (2.32) 

J lIM-pY J 

+ usual helicity-2 contributions from GR . 

Here we have used the fact that the field <fi can be expressed in terms of the Feynman 
propagator 

(j){x) = — f d 4 x'G F (x,x')8T(x'), (2.33) 
V6Mpi J 

where we have defined the Feynman propagator 

UG F = i8 4 (x - x) , (2.34) 

and where the modified d'Alembertian operator □ is defined in (2.25). 

As usual the Feynman propagator can be expressed in terms of the Wightman 
functions 

G F {x, x') = 6(t - t')W + (x, x') + 6{t' - t)W~{x, x') , (2.35) 

where 

poo 

W + (x,x') = J2 / du uMuU^YU^^Ylje^^e-^-^, (2.36) 

and the mode functions expanded in spherical harmonics U£ mul (r, Q, t) = u^{r)Y^ m {Q)e~ wt 
are the eigenfunctions of the mode equation Du£ mul = and form a complete set. 
The time averaged power is 

. jp . poo 

P = -(-^)= / dww/(a;), (2.37) 



\ dt I j 

where /(a;) is determined from the imaginary part of the effective action 

2 f°° 
lmS cS = / duf{u) . (2.38) 



Tp Jo 
We define the moments 



M tmn = ±-jJdtjd?x u ln {r)Y lm {6, ip)e- mnpt 5T(x, t) , (2.39) 

where we use the notation ug n = ug )UJ =nn- Taking the Fourier transform 

00 

M em = M lmn e inQp \ (2.40) 
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we have 

f(u) = r dt f dt ' Re U^MU^MlJt')) , (2-41) 

where we have used the facts that the imaginary part of ue m integrates to zero and 
that Ylm=-e Yi m (9, <p)Ytm(Q' ■, if') is real. Then the period-averaged power emission is 

oo 

P1 n=0 im 

Notice that since we only integrate over positive frequencies in (2.37), we only need to 
sum over positive harmonics n > 0. 



3 Power Emitted in the Quartic Galileon 

First we consider the case that there is a single strong coupling scale, A 3 = A 4 = 
A 5 = A, so r*^ = 7* 4 = r*. We compute the power emitted in the Galileon using 
the formalism developed in the last section. This requires us to derive the properly 
normalized mode functions ii£ n . 

3.1 Mode functions 

The field fluctuations can be expanded in terms of the mode functions 



t)= J2 a imnUi n (r)e- inQpt Y im (n), (3.1) 



n=— oo £m 



where we have used the fact that for periodic systems we need only sum over a discrete 
set of harmonics n instead of integrating over a continuum of frequencies. 

The radial mode functions ut> n are solutions to the homogeneous equation 

Uu ln e- mnpt Y lm {9) = 0, (3.2) 

subject to the normalization defined by equations (2.34) and (2.36), which is valid as 
long as the field reaches the oscillating WKB regime within the strong coupling region 
(i.e. as long as R+ > Q" 1 ). In practice however we compute this normalization by 
matching with the free Minkowski spacetime normalization at spatial infinity, r ^> r* 

1 cos(nQ P r + P) 
hm u in (r) = —== , (3.3) 

where the phase P is irrelevant. 
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3.1.1 Strong Coupling Regime 

In the strong coupling region where £4 dominates, rCr^ the mode functions which 
satisfy the correct boundary conditions at the origin are (see Appendix A for a detailed 
discussion of the choice of boundary condition) 

u 0n (r) = u 0n cosuj 0n r (£ = 0) (3.4) 
Ui n {r) = ue n sinu en r {£ > 0) , (3.5) 



where 



c^l(nfip) 2 -^^^, 



and where the normalization u is fixed for each mode by matching onto the correctly 
normalized mode at infinity according to (3.3). 

Matching the WKB solution with the strong coupling solution at r = r* we find 



1 



£ -C nL 



, e n , i > n£ CTit 



where £ cr i t is defined by 



-crit 



Q P r* . (3.8) 



Modes with imaginary ue n are exponentially suppressed, so when computing the power 
we only need to sum over modes with real ue n . With this in mind we use the following 
approximation in what follows, 

u in = d(n£ clit -£) __ , (3.9) 
V7rm lp 

where we define the step function as 9{x) = for x < and 1 otherwise. 
3.1.2 General Form of the Power 

As usual, by conservation of energy the monopole does not radiate in the non-relativistic 
limit [i.e. at leading order in the velocity expansion. See Ref. [35, 36] for the behaviour 
at next order). Similarly, by momentum conservation the dipole does not radiate at 
leading order in the velocity expansion 5 . The power for the higher order multipoles is 
given by 



[ P dt [ d 3 x e~ innpt sm(nQ P r)Y em (9,^)5T( y x,t) 
Tp Jo J 



2 



5 Normally one associates the dipole moment with the moment linear in the velocity v and by 
momentum conservation this moment does not radiate. In our case whilst it is true that the dipole 
moment defined as £ = 1 does not radiate, the higher order multipoles also have only one power of i>, 
and they do radiate. The spherical harmonics with £ > 1 allow for time dependence. 
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where the expression for the perturbed source 5T(x, t) is given in (2.18) and is propor- 
tional to 5^(f—ri(t)), and the two objects follow the standard Keplerian orbits fj(£), 
given in spherical coordinates by 

ri , 2(() = ^S^h. < 3 - n > 

M tot 1 + e cos 1 l pt 

71 

with 6* lj2 (i) = — and ^(f) = fipt + , (3-12) 

where e is the eccentricity and f is the semi-major axis of the orbit. For simplicity 
(and without loss of generality), we choose the plane of the orbits to be localized in 
the plane 9 = it/ 2. The remaining angle ip is then determined knowing that the two 
objects orbit with frequency Qp and are always diametrically opposed on their orbital 
path. 

Defining the reduced mass 

MiM 2 

■^iP' (3 - 13) 
then to leading order in Qpf the power emitted is 



M 2 n 2 P f 2 

2 ~,2 

PI * n=0 im 
1 + ( — I)"" jj. „-i(n-m)n P t 1 — e 2 



« = ^^EE^--**)" 2 ^ (i,o) (3.i4) 

^ e -i(n-m)fl P t _ 







e cos f2 pt 



Note that in Minkowski spacetime radiation problems, the mode functions are je(oor) rs 
(wr) at small distances, and so higher order multipoles are suppressed by more powers 
of Qpf ~ v. Here however there is no additional velocity suppression for higher order 
multipoles (for as £ < n£ cr i t ). 

The integral can be evaluated as 

1 r Tp p -i{n-m)Clpt o_ / n n-m 

^ / dt ri — o7 = (- 1 ) n " m ^=^ ; — 5=i > ( 3 - 15 ) 

T P J 1 + e cos Qpi V 1 - e 2 \ 1 + V 1 - e 2 / 

for n — m > 0, so the power is 

(PH^Vm'n.*., ( 3, 6 ) 



M 



PI 



with 

«4rit ^ 2 / \ 2(ra-m) 



n=0 £=0 m=-£ 



2 fmn\ 

y/T^J V 2 / 



cos 



Note that the £ = 1 mode does not radiate (since ^1^(^/2,0) = 0), as expected 
from momentum conservation. So the first multipole that has nonzero radiation is the 
quadrupole t = 2. 
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3.1.3 Quadrupole 

Let us now compare the power emitted in the quadrupole to the power from GR and 
from the cubic Galileon (in the case where the quartic and quintic interactions are 
absent, but keeping the same strong coupling scale A). The power emitted by the 
quadrupole in the n — 1 harmonic is 

(P)W - ( M V {Qpf)2 Q 2 (3 17) 

Comparing this result with that of the cubic Galileon presented in Ref. [35], 

(^P) Cubic Galileon ~ J (OpT*)^ J2 P ' ^ 18 ) 

We see that the relevant Vainshtein screening is (fipr*) -2 compared to the Vainshtein 
screening appropriate for the force between the pulsars (f/r*) 2 . This is exactly anal- 
ogous to what happens in the cubic Galileon, where the Vainshtein screening is less 
effective than in the static case. We also see that the full Galileon is enhanced by a 
factor (Qpf)^ 1 ~ v" 1 relative to the cubic Galileon because the quadrupole is sourced 
by the monopole moment. 

So far, we have only considered the contribution from the first harmonic. The 
Galileon radiation includes radiation from all harmonics, each of which contributes 
equally to the moment. 

3.1.4 Summing over all multipoles and harmonics 

If one were allowed to sum over all harmonics (till n — > oo), the power emitted would 
formally diverge. In order to gain a better understand of this divergence we truncate 
the sum over the harmonics at the cutoff n max in (3.17), and denote as ^(fimax) this 
regularized sum. Furthermore, for simplicity we focus on the case with no eccentricity, 
i.e. e = 0, so that the regularized sum simplifies to 

So(rw) = £ E E ^ cos2 (t) • (3 - 19) 

n=0 £=n m=—£ 

Since Qpr* ~ 10 6 ^> 1 for realistic pulsars, most of the terms in the sum have i ^> n. 
We can then use the approximation valid for £ ^> n, 



) « ^cos((£ + n)|) . (3.20) 



The details of the calculation can be found in Appendix C, and the final result is 
the power depends quartically on the cutoff n max . 
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To get a sense of this result we apply it to a pulsar system consisting of two 
solar mass objects orbiting with a period of 2n/Qp = 8 hours with a semi-major axis 
f = 10 9 m and with no eccentricity. This choice of parameters is close to those of the 
Hulse- Taylor pulsar [41]. 

The GR result for this system is given by the Peters-Mathews formula (assuming 
zero eccentricity) 

2 5 M\MlM tot 



Peters— Mathews 



5 r 5 

Comparing this to the naive Galileon result (3.18) we find 



(3.22) 



Ppull Galileon _ g x 1Q _ 4 ^ (3 33) 

-^Peters— Mathews 

Even for n max ~ 1 this is a large amount of power compared to the cubic Galileon case 
in Eq. 3.18. However if we trust this calculation then there is no natural cutoff in n 
until n ~ (Vtpf)~ l when the assumption we made that wf < 1 breaks down. If we 
take this cutoff and use Hulse Taylor parameters we find that ^tp^~ ~ 10 9 . In what 
follows we will interpret this result as a breakdown in perturbation theory. 

3.2 Validity of Perturbation Theory 

The divergent power suggests that our calculation was too naive. Since we have been 
using linearized perturbation theory, the natural thing to check is whether the fluctu- 
ations themselves become nonlinear. Indeed, we might expect perturbation theory to 
break down on physical grounds. We have used perturbation theory around a spheri- 
cally symmetric source, but have found important contributions from arbitrarily high 
multipoles £. Since higher values of £ are more sensitive to what happens over small 
angles, we expect that our choice of background should become worse for large I. 

In this section we check the validity of perturbation theory around the spherically 
symmetric background, by explicitly constructing the first and second order pertur- 
bations. Since we want to compare the physical values of the fields, we use here the 
retarded propagator. 

Based on the discussion above we cutoff the sum over I in the propagator. If one 
can trust the perturbation series at all, it can only be trusted at low £. Physically this 
is because there is some uncertainty associated with the angular position of the two 
objects. Trusting perturbation theory to arbitrarily high £ and performing the sum 
over all £ would imply that the positions are known with arbitrary accuracy. High £ 
modes are also more sensitive to the non spherical nature of the source, and at high 
enough £ we do expect the assumption of a spherical background to break down. Thus 
we introduce a cutoff L on the sum over £. We take the cutoff L < £ crit , since £ crit is 
itself very large for realistic pulsar systems. 

The analysis is performed in the WKB regime where derivatives acting on 
can be expanded in powers of fip . As we shall see below, the field diverges at certain 
isolated points and we compare the field values at these points because they will give 
the largest values of /(j)^. Finding this ratio to be bigger than one at a single point 
is sufficient to show that perturbation theory is breaking down. 
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The exact equations of motion for the Galileon ir(x) is 
□tt + j- 3 ((Dtt) 2 - {dM*) 

+i [(Dtt) 3 - SDtt [d^f + 2 (S^tt) 8 ] = . (3.24) 

Performing a background/perturbation split for the source, T = To + <5T, we will be 

interested in the second order perturbations, for the field it = n(r) + a/2/30^ + 

2/30( 2 ) + • • ■ , with T ~ 7r(r) and 5T ~ as already computed in the previous 

section, (in particular in the strong coupling region the background configuration for 

7r(r) is given in (2.17)). The second order fluctuation is sourced by nonlinearities 

in 0^. To check the validity of perturbation theory, we will compare the magnitude 

of (2) to the magnitude of <p^\ The perturbative expansion is under control only if 
0(i) > 0(2) f or a u r> 

Unlike in the previous section where we were interested in the power emitted, 
we here focus on the physical values of the fields, and thus construct the retarded 
propagator, 

G R (x,x') = -i8(t - t')(0\[(f)(x),(j)(x')]\0) (3.25) 
= -9{t - t') J du sin u(t - t') (3.26) 

L I 

E uar)uUr')YU8,v)YL(9',<p'). 



L i 
1=0 m= 



In the WKB regime, the first order field fluctuation ft 1 * is then given by 

</>«(*) = - J ^xG R {x,x')^l (3.27) 
2nMf j^^^-mi^) (3.28) 



x Ye m (9, 0) Yi m y— , 0J sin mflpr sm(mQpt — mcp) cos 2 j ? 

where the sum over the multipoles has been truncated at the cutoff L, see Appendix D 
for details of the above calculation. 

The fluctuation reaches its maximal value on a radial light cone when r = 
t = (2k + l)n/2Qp with k e Z, with 9 = ir/2 and <p = 0, (note that m is then forced 
to be even). Calling this set of parameters x max (note that x max is not a unique set of 
parameters, it is just any choice on the radial light cone that satisfies these conditions). 
If one were to push L — > oo, the sum for <p^ would diverge exactly at x max . However, 
in terms of the finite cutoff L, the value of the field at x max is given by 

m , ,1 2nMf L 2 4 M f „ , 

(1) x max , L) « — -= - = = — — -L 2 , 3.29 

V V max ' ; 4.5 y/6M Pl r 2 vr 2 9iry/E M Pl r* V ; 
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where the factor of 1/4.5 is an approximation 6 . 

The crucial observation is that <fP-> does not fall off with increasing r. This 
is ultimately tied into the fact that the effective metric for the fluctuations is one 
dimensional for r<r t . The result of this is that (p^ does not see a compact source of 
size f or even Up 1 as might be expected, but rather one that extends out to r*. This 
pumps a large amount of energy into the second order perturbation, and this is what 
ultimately makes perturbation theory break down. 

The perturbed equation for (j)^ is 



°</> (2) = in 



A 6 



(3.30) 



5u^7r o ^^0 (1) D0 (1) + 2d,d u 7T d u d x <f) m d x d' 1 ^ 1 



1 

A 3 



(□<^>) 2 -(^«) 



where □ is given by (2.25). 

The source for the second order fluctuations is in principle complicated, but we 
simplify its expression by working in the WKB regime r > Up 1 , and focusing on the 
leading terms in flpr. The first order field fluctuation then takes the form 



(£, t) ~ A(r)B(9, if) cos(nQ P t + P t ) cos(nQ P r + P r 



(3.31) 



where A(r) is a slowly- varying function of r (which varies over distances much bigger 
than Up 1 ) and P^ r are irrelevant phases. For I < £ crit = fipr*, we have ~ 
d 2 <p^ ~ J2 n ( n2 ^p + Z T £ + ' ' ' ) < P^ ^ ^Vn^ 1 ^, so we can ignore the angular deriva- 
tives and focus only on the leading order contribution from the radial and time deriva- 
tives. Recalling that d 2 n ~ M 1 / 3 A 2 /r, the different contributions sourcing 111 
(3.30) are then of the form 



A6(^o)(<9 2 < 



(ih2 



^M^A 2 
A 6 r 

-- (r 

A 3 r ^ 



2 /2o4 , ?7-2ri ^P , 

n n U P H h 

r 



(0 (1 



(3.32) 



and ^(9 2 0W) 2 




2 /2r>4 

n n U P 



(3.33) 



The 



-^^f2p</> 2 contribution in (3.32) should clearly be the dominant one, however a 
straightforward calculation show that it actually vanishes exactly, and the contribution 
from (3.32) is thus of the same order as that of (3.33). 

= T^^d 2 from the cross term in 

A J UJT Z ~ 



The next order corrections are the 

A3 



1 r*ar 
A3' 



3.32 (that arises from £4) and the T3</> 2 contribution coming from 3.33 (that arises 



6 Naively one would guess this factor is 1 /4 because we are taking roughly half of the m terms and 
half of the I terms. Numerically it seems that 4.5 for this factor is a better fit. 
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from £ 3 ). We see that at small r, the £4 contribution is dominant (as expected), but 
becomes subdominant when r > (r^/co) 1 ^ 2 = -^(VLpr*) l l 2 VL~ I } = -^(Qpr i <)~ 1 /' 2 r+, i.e. 
still within the strong coupling regime but already in the WKB region. This is true 
for all n. Since we are interested in evaluating the source at r = r*, the contribution 
from £3 is the most significant. This is explained in more depth in Appendix D. 

We can solve the equation for by using a WKB-like ansatz. In the limit 
L —7- 00, we would find again a diverging expression on the radial light cone, however 
keeping a fixed cutoff L we find 

0( 2 )(x max ,L) = ^VAf ^ Mf \\l— A . (3.34) 

We can then explicitly compare <fP-' and (j)^ 

L) 3V6^ M 1 fip r4 
<j>W(r*;x max ,L) ~ ~m~Wi V WJ^ 1 ' J 

= JL<[IJ± v n pn L 4 (3.36) 
968tt V2M tot v ; 

^O.lxt^L 4 , (3.37) 

with the orbital velocity v given by v = Qpf, where we assumed A4. ~ M tot , in the last 
expression, i. e. two bodies of comparable masses. 

We find that in principle there are systems where perturbation is valid over some 
range of L so long as the system is sufficiently light or slow (small vQpr+), i.e. as long 
as 

(M V 1/4 
L < i^vVLpn ■ (3.38) 



tot- 



However for realistic binary pulsar systems, vQpr* ~ 10 3 , so perturbation theory 
breaks down for all I. 

3.3 Effect of a Multipole Cutoff on the Power Emitted 

As mentioned previously, one can only trust perturbation theory as long as the second 
order field fluctuations are small relative to the first order one (in principle this is not a 
sufficient condition, but it already gives a good handle on the behaviour of perturbation 
theory in a given setup). We thus consider a cutoff L which small enough and replace 
5 (n max ) in (3.19) with Sq{L), that is to say we replace the cutoff in n with a cutoff in 
£ with L < £ crit . Recomputing the expression for the power with this cutoff, we get 

1 00 L £ 2 r4 

5«W = ^EE E n ^ (I' °) <W - ^ , (3-39) 

n=0 £=1 m=-£ 

where the approximation is obtained numerically 7 . The factor of m? weights the high 
m modes (where m « i more heavily and so this approximation breaks down). In 



7 If we had used the approximation t ^ m we could have done the sum exactly and would have 
found L 2 /487r 2 . 
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terms of the power emitted this implies 



(P) = ( 



3 M 2 , \r, 

4 M 2 v 2 
27M 2 Y (fi P r + 
4 121 M 2 



ftp! 4 (3.41) 



n 2 P . (3.42) 



Note that this result is more Vainshtein suppressed than any one mode (Ptot ~ r * 3 
instead of P mo dc ~ r * 2 )- 



4 Hierarchy of Masses 

We now investigate setups where perturbations theories is under control. One of the 
requirements for that is given in Eq. (3.38). In particular we see that in the limit where 
one mass is much bigger than the other, then 

M m 2 «m^ M. <<; 1 



Mot M 1 

and so we expect perturbation theory to be under control in that case. Taking the 
upper bound for L, L = (-^-t>f2pr*) -1 / 4 the power radiated in that case is then given 
by 

(P) « 100 x ^ j^L n 2 P . (4.2) 



Mtot (fi P r, 



As a fiducial example of such a system, let us consider the Earth/Moon system and 
use the power emission to estimate the rate of change of the orbit to check that we 
get a physically reasonable result. In this case using M Earth = 5.97 x 10 24 kg, M Mo on = 
7.35 x 10 12 kg, v/c = 3.42 x 10~ 6 , and Qpr* = 648 we find perturbation theory works 
with a cutoff L = 73. Since the eccentricity of the earth-moon system is e = 0.05 the 
approximation of zero eccentricity we have made in the calculation of the power is a 
good one. We find a radiated power 

(P>Moon « 10- 112 M|, , (4.3) 

which implies a rate of change of the orbit of 

f d-E/NR 

r = 

-Enr dt 

where the system's non-relativistic energy is 



(4.4) 



i m x m 2 ( n 2 \ 1/3 . . 

^nr = 9/ , I TT-TT- M ?i ■ (4-5) 



87T 2 /3 Mpi V^totMp! 
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Comparing this to the GR result for the power emission using the Peters-Mathews 
formula (assuming e = 0), 



Peters-Mathews 



32G 5 M|M2 oan (M e + M ] 



Moony 



5 c 5 r 5 



(4.6) 



we find 



^-~1(T 33 , (4.7) 
rcR 

which is utterly negligible ! We see that in that setup, the Vainshtein mechanism in 
the quartic Galileon is very active and prevents the scalar field from radiating almost 
any energy from the system. 



5 Hierarchy between Two Strong Coupling Scales 

A second situation worth mentioning where perturbation theory can remain under 
control in a binary system in the presence of quartic or quintic Galileon interactions is 
when these interactions do not dominate straight away but only far within the strong 
coupling regime. To be more explicit we consider in what follows a hierarchy between 
the different strong coupling scales A3 <C A4 (or equivalently r*^ <C f*,3) 8 , so that as 
one probes shorter distances, the cubic Galileon starts dominating first and only very 
deep in the strong coupling region does the quartic Galileon take over. 

To explore this situation, we construct the mode function in stages. We first 
approximate as being equal to its leading order behavior in each region. We start 
with the £4 region, which extends from the origin out to r* >4 . We first need to determine 
the correct boundary condition for the mode at the origin. It is not necessarily correct 
to take the mode that is smooth at the origin, because the equation of motion is singular 
there. So long as we are are in the region r <C ^p 1 , the field can be approximated 
as a power law, and we can extend the solution beyond the 'crossover radii' 7^3 and 
r^4 by matching <fi and its first derivative at these radii. However once we reach the 
regime r > flp 1 , the field begins to oscillate and we can no longer use this matching 
procedure. Instead we can use the WKB approximation to extend the solution out to 
infinity, where we can then fix the normalization. 

It is clear that the details of this procedure depend on which strong coupling 
region the scale Hp 1 falls into. To get a sense of typical scales, take A3 = (1000km) -1 
and consider two solar mass binary pulsars with orbital period of 8 hours (which 
approximately describes the Hulse- Taylor pulsar), then Qpr* 3 ~ 10 6 . So the case where 
fip 1 <C r^ 3 is of most physical interest. In what follows we will further assume that 
r* )4 <C dp 1 <C r* j3 , so that the crossing between slowly- varying and the WKB region 
occurs in region where the cubic Galileon dominates. In the case where 4 3> fip 1 , 
then the result is almost insensitive to the presence of the cubic Galileon. 

8 Of course we still assume r*,4 3> f, otherwise the quartic Galileon would not start dominating 
before one starts being strongly sensitive to the internal structure of extended object itself. 
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5.1 Mode functions 

The modes in the £ 4 region are again ui n = ui n sin(uj£ n r), for £ > (see eq. (3.5)). We 
then match this to modes in the £3 region at r = 7**3 where the modes are given by 
[35] 

f r \ l,A V3 / r \ V 4 ^3 
utn = a ln ( J Vl (—ur) + b ln ( Y Ue (—ujr) , (5.1) 

V r *,3/ ^ V r *,3/ * 

with vi = (2£ + l)/4. Matching U£ mn and its first derivative at r = r*,4 <C fip 1 , we find 
for £ > and to leading order in lot* 4 



a fa = « to l±ir(i* + 1) V4 f ^ 1 (nfipn^) 1 -' (5.2 

6 fc = 0((nfi P r M ) 1+ ^). (5.3) 

Now one can finally match this solution to the WKB regime in the limit flpr ^> 1 
to fix ug n . Since b^ n ~ 0, one can simply use the results derived in [35] and set 
a ln = (9n/128rl 3 ) l /\ so 

V4 r/o\ vt ^o_^ ^^-1^/4 



1 + 2£ / 9tt \ 1/4 / v 7 ^ \ (nQpn^- 1 



U£n 2 + £ [l2s) [a) ' T(V+ 1) r£ 4 (5 ' 4) 

/r ,\ 1/4 1 

= PtinQpr^) 1 "- 1 -— , (5.5) 

\ r *,3/ 

where = (^) (rll) 1 ^ 4 + * s a dimensionless prefactor. Note that 
is order 1 for £ = but falls off rapidly with large £. 

5.2 Power 

Using these correctly normalized mode we find for a zero eccentricity source that 

o / 00 

<*> = ^rrSS" 2 " 1 ^^ 2 ^ 1 ' (5 - 6) 



f^l P dt e ~ Mpt J d * x Mutnr)Y lm (6,<p)6T(x,t) 



For £ <C n£ CV it we have ~ nQp. Since Opr < 1 we take the leading order behavior 
sin(nf2pr) « nQpr, and find 

3M P1 (fi P r +i3 ) n=0 lm * 



Consider the sum 

00 

n^fj^Upr^y^-^YeJ 71 

n=0 



S(£) = J2 n^ +1 $(n P r* 4 ) 2 ^Y en (l, 0) 2 • (5.J 
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Here the small parameter fipr*^ plays the role of the velocity in the cubic Galileon 
by suppressing the higher order multipoles. At large £, using Stirling's approximation 
T(z) ~ \/~2ixz z ~ x l 2 e~ z and approximating the sum over n as an integral, we find the 
scaling behavior with I 



eV3 



s{*)~ l-^-JW) r 1 ' 2 . (5.9) 

Thus for flpr^ < \/3e/2 « 2.3 the summand becomes exponentially suppressed at 
large t. At very large Qpr+^ the sum is dominated by the low multipole. Numerically 
we find that for Qpr*^ = 1CT 2 that S(4)/S(2) « 1CT 3 , and perturbation theory is then 
well under control as we recover a suppression at higher multipoles and moments. 

5.3 Quadrupole Radiation 

Focusing the power emitted by the Galileon quadrupole i = 2, we have 

oo 2 2 2 

\fi E E nV2y ^ (i> °) «w = ~S x 2?/2y22 (I' °) w °- 1 ' (5 - 10) 

n=0 m=— 2 

and so the power in the quadrupole is 

This result is almost identical to the cubic Galileon result in Eq. 3.18, but with one 
power of velocity replaced with Qpr*^. The power is velocity enhanced relative to the 
pure £ 3 case because the £4 modes pick out the dipole moment of the source. 



6 Discussion 



We find that the behavior of perturbations around spherical, time dependent back- 
grounds depends strongly on the presence of the fourth Galileon interaction. Unlike 
the case of the cubic Galileon considered in [35, 36] we find that studying perturbations 
around a naive static, spherically symmetric background is insufficient for computing 
the power emitted by a binary pulsar system. The key difference is the effective one 
dimensional metric seen by the fluctuations 

Z^dx^dx" oc -dt 2 + dr 2 + r 2 4 dfi 2 . (6.1) 

As a result of this one dimensional metric, the naive perturbation theory predicts that 
all multipole modes with fixed n (defined by u = nVtp) contribute equally to the 
power until I ~ nQpr*. Since for typical pulsar systems Qpr+ ~ 10 6 this means a huge 
number of multipoles radiates with comparable strength. For simplicity we focused on 
systems with orbital eccentricity e = 0, but the divergence would be present even if 
we included the effects of eccentricity. 
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The resolution is that the 'spherical background plus non-spherical perturbations' 
approximation breaks down. Physically we expect that this occurs because modes 
of arbitrarily high I contribute to the power emitted, and higher i modes are more 
sensitive to the lack of spherical symmetry in the system. This was checked by explicitly 
constructing the first and second order physical solutions to the equations of motion 
using the retarded propagator with a cutoff L in the sum over multipoles, and taking 
their ratio at a point where their ratio was maximized. We found that whilst there is 
a range of parameters where perturbation theory could be trusted up to some L, for 
realistic pulsar systems we are forced to take L < 1 so naive perturbation theory is 
never valid. 

We expect that this result implies that there is an additional Vainshtein screening 
on top of the usual static and spherically symmetric screening. The breakdown of 
perturbation theory indicates that the perturbations themselves are nonlinear, and so 
contribute to their own Vainshtein suppression, in addition to the normal Vainshtein 
suppression from the background. As preliminary evidence in this direction, we note 
that the expression for the power with the cutoff in L is more Vainshtein suppressed 
than the power calculated from any one mode. 

We also studied two regimes where perturbation theory can be recovered, to check 
that the theory makes sensible (small power) predictions in these cases. The first 
situation is when there is a hierarchy of masses between the two bodies in the system 
(such as the Earth-Moon system). As a second example, we consider a hierarchy 
between the strong coupling scales of the cubic and the quartic Galileon. The WKB 
oscillating behaviour then starts in the strong coupling region where the cubic Galileon 
dominates over the quartic one. We then find a suppression in the higher multipoles 
as expected from previous results. Nevertheless, as long as the interactions from the 
quartic Galileon are important on scales comparable to the size of the system f, the 
system radiates as if it had a size r*^ rather than f. 

We expect this result to be a generic feature of time dependent systems exhibiting 
the Vainshtein mechanism when the quartic Galileon is included. Thus the intuition 
gained from studying the static, spherically symmetric case may not directly apply to 
more complicated time-dependent (or less symmetric) systems. This two body system 
is surely the simplest generalization of the one body, static, spherically symmetric case 
and we already see at this level that a more detailed understanding is required. 

Future work should involve going beyond this naive approximation to get a better 
analytic handle on the power emitted. We expect that one can find a background that 
takes into account the time dependent evolution. 
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A Normalization of the Mode Functions 

In this section we check explicitly check the mode normalization by expanding around 
a regularized background. Consider a smooth background source no longer localized 
at the origin, but rather spread over a radius e, 

ir = ^-r)^«, (A.l) 
where e <C r*^. The background for r < e is given by 

E vr'(r) A 2 /M/M P1 \ 1/3 



Af / M/Mpi V 
e V 24tt J 



(A.2) 



The equation of motion for the fluctuations valid for r < e is then 



6 /M/M P1 \ 2/3 1 



A 2 . V 24vr 7 e 2 _ 

2/3 



-3^+la r (r 2 9 r 0) + lv^ 



0, (A.3) 



so that Z rr = Zqq = \Z tt = -£i ^ M 2^ P ' j ^- Now there is no singularity at r = 0, so 
modes inside r < e satisfy 

u fmn( r ) = utmnjeiVSur) ft u imn (V3ur) e . (A.4) 

Outside the source, for r > £ the modes are given by 

u imn = A cos Logr + B sin u^r ft A + £?a;r , (A. 5) 

using the approximations £ <C flpr^^ and r -C r +j4 . 

Now the equations of motion are smooth at r = e, and so the appropriate matching 
conditions are M^ n (e) = < mn (£) and 9 r w^ mn (e) = d r uf mn (e). Solving these for A and 
B and taking the ratio, we find 

B = —" £ ' (A ' 6) 

so in the limit e — > the term going as B dominates over the constant going as A (this 
is also valid for the dipole £ = 1). 

For the monopole (£ = 0) on the other hand the contributions are different and in 
that case the same analysis shows that the constant term A dominates. This is reason 
why the mode normalization differs for I = and £ > as can be seen in eqs. (3.4) 
and (3.5). 



B Small Strong Coupling Scale for the Quintic Galileon in £3 
Region 

Here we find that Z rr and Zqq are dominated by the pieces that are independent of 
A 4 and A 5 , and so these reproduce the equations of motion for fluctuations from the 
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cubic Galileon. There is a contribution from £5 to K t however: 

32 /rv 3 \9/ 2 MsV 72 , e A /r*, 3 \ 3 / 2 



+ higher corrections . 

The higher corrections piece refers to terms that are higher order in r/r*^ and also 
terms that are suppressed by powers of A3/A4 (which must be a small ratio for this 
region to even exist) which do not dominate for r > 7*4. 

The first term here is comparable to the second term for r ~ (A^/ A^A^r*^. Thus 
if we take A 5 > (A 4 /A 3 ) 1 / 3 A 4 then we can ignore this term in the £3 region. 

If we do take a small A5 then we find that the equation of motion has the form 
of Laplace's equation in 3D for r < (A4/A3A5). 




1 /r*- A 3/2/' 2 . , , I 16 - /A 3 



9/2 



' (B.2) 

This is similar to what occurs in the £4 region, where time and angular derivatives 
appear without a relative factor of 1/r 2 . However a detailed analysis of this situation 
is beyond the scope of this work. 

C Sum over the Multipoles and Moments 

We start with the sum in (3.19) 

Sb(wmax) = ^2 Y2 X] n2 ^ m (f '°) cos2 Sn ' m ' ( C ' 1 ' ) 

n=0 £=1 m=-< 

Now focusing on multipoles with I ^> n (for a realistic pulsar system, Qpr* ~ 10 6 ^> 1), 

F^(|,o) «icos((£ + m)|) . (C.2) 

Note that F to (f,0) 2 cos 2 (^f) « ^ cos 2 ((£ + m)f ) cos 2 (^f) = ^ cos 2 (f ) cos 2 (^f). 
Then the sum becomes 

So(rWx) « ^2 X] ™ 2 cos2 ( y ) cos2 ("Y") 5n ' w ' ( C - 3 ) 

Now we can figure out the value of the two inner sums as a function of n. Doing 
these sums amounts to counting the number of nonzero terms, because the summand 
is either 1 or 0. Explicitly 

n^crit I / ft \ 

(C.4) 



S (n) = n 2 cos2 ( y ) cos2 [~yj 5n - 

l=\ m=-£ ^ ' 
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Then we arrive at a final approximation for the sum 



So(rw) - g E - 3 cos2 (f ) « 7^<a* • (C6) 

n=0 

D First and Second order field fluctuations 

D.l Retarded Propagator 

We normalize the retarded propagator so that 

nG R (x,x') = 5\x-x'), (D.l) 

where □ is given in (2.25) (strictly speaking this normalization procedure is only ac- 
ceptable as long as the modes are continued all the way to the oscillating WKB regime). 
In this procedure the mode functions are real, U£ W {r) = u* iuJ {r) and so we have (using 
the fact that Y£ m Y e * m is real) 

/oo I 
dwsinw(t - t') ^ u £Ld (r)u eul (r')Yf> m (6, ip)Y e * m (6\ <p') . 

£=0 m=-l 

(D.2) 

Then given a source J(x) (D0 = J), the field is given by 



<f)(x) = j d 4 x'G R (x,x')J(x') (D.3) 

= - J dV J du9(t - t') shiest - t') (D.4) 

oo I 

Y Y U£ u] {r) U £ LU (r')Y£ m {e,if)Y; m {e\if')J(x') . 



oo t 
£=0 m=-£ 



Now define 

■WO = J d 3 x'uUr')Y; m (8', V ')J(x',t'), (D.5) 
and since the source is periodic we can write 

•W*)= Y fn(u,£,m)e mnpt . (D.6) 

n=— oo 

Then the expression for the field is given by 

- Jdt' J dtu9(t-t')smcu(t-t') Y Y u ^( r ) Y trn(OA)fn(co,e,m)e mnpt ' 

n=—oo im 



[X 



(D.7) 
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Performing the integral over t' yields 

oo „ 

<t>( x ) = E E / duj e mnpt f n (ooJ,m)u eu] (r)Y em (0^) (D.8) 



-oo Im 

-— (S(u — nVtp) — 5(u> + nVtp)) 



2 v ~ v ~ " v ~ ' '"""" yy ' u 2 -n 2 Q 2 P 

= M*) + Mx), (D.9) 

where (pp(x) is the piece with the delta functions and 0s (x) is the piece without the 
delta functions. As we will see, <pp is the radiating piece of the solution and dominate 
in the WKB regime while <fis is essentially the static piece and is negligible in the WKB 
regime and does not contribute to the radiated power. 

We first focus on simplifying (pp. The u integrals are trivial because of the delta 
functions. Carrying through the simplification we find 

(j> R (x) = 7r^5> fa (r)y £m (^) / — / tfx'smnn P {t-t')YZ m {e\v>in{r')J{J,t')- 

n=0 Im J P J 

(D.10) 

Note that for Qp = 0, (ftp = 0. This part of the field is only present for time dependent 
sources, and it is responsible for radiation. We could have seen this above too, we had 
5{u) — nVtp) — 5{uj + nVLp) = for Qp = 0. 
As for the second part of the field (ps, 



oo „ 

E E / ^e innpt u u Y lm {0, <f>)f n (u, A m) ^ _ " (D.ll) 



n=— oo im 

oo 



dt' f 3 , f , UJ 



= 2 y ^ y dV y ^ ^- (napY Y ^ Y ^^ < d - i2 > 

n=0 £m ^ ^ J J y ^ ' 

^u u t(r)u u i(r') cosnQp(t — t')J(x", t') . 

Note that this contribution does not vanish for Vtp = 0, so it includes the contribution 
from the static propagator. 

This term is suppressed by compared to (f)p as can be seen by the following 

argument. We do not expect the n = mode to contribute to the power because the 
static piece is not oscillating so derivatives acting on it are suppressed compared to 
the oscillating parts of the field. Note that for large n the integral over u is basically 
because u U £(r)u U £(r') is even in u (at least for all of the kinds of modes we have 
considered) and u/ (u> 2 — Uq) is odd in u. For large n, the pole is far enough away from 
that the bounds on the integral can be approximated as , which gives because 
the integrand is odd. 

So since n = and large n modes are small compared to the radiation piece, we 
expect that this term should be small in the radiation zone. We can explicitly see how 
this works for the £4 modes by using the integrals 

/ \ . /1 \ ( cos(ap) sm(bp) if a > b > 

sm ax sin (cm )r , ; / r ' .. 

' ±sin(2ap) if a = b > 

sin(ap) cos(fep) if b > a > 



00 



dx v ' : ; =--{\ sin(2ap) if a = b > . (D.13) 
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If we use the modes appropriate for £ 4 for < t <C £ CT it, that is u n = ^-^-^ sin(72f2pr), 
then the integral over u in D.ll is of this form with x = u, a = r, b = r', p = nflp. 
Then we can see explicitly that for fl P r 3> 1 that (fis is suppressed with respect to (f)R 
by 1/nQp. 

D.2 Field in the WKB regime 

Now we explicitly construct in the WKB regime. We consider only the radiating 
part, (and drop the subscript). The source for 0^ is J^\x) = — -5T(x). 
Here we plug in the modes appropriate for £ 4 . 

oo 

bV>(x) = 7rJ2J2 u ^ Y t™( e ^) (D.14) 

o em 

^ J d 3 x'smnn P (t-t')Y; m (6',<p')u in (r')J@,t') 
Z 2 EE M J T,(l + (-l)nsinnfipry, m (^^)y, m (^0)(D.15) 



n=0 Im 

X 



X 



d/ 7 

— 9(£ - n£ CTit ) sin nfip(t - t')e~ mnpt f e (t') , 



where 



ri,2 = 1 r 7 J -r. (D.16) 



Mtot 

We define the reduced mass 

x - M X M 2 M 2 M 1 , , _ , 

E^-^-+^ = ^-, (D.17) 

where M. is the reduced mass. Considering a source with no eccentricity, we have 

kit') = i. 

Then the integral over t' simplifies to 

y ^sinr^-^"™^^ . (D.18) 

We can simplify the resulting expression for <p (remember n > and use the Shortley 
phase rule Yt- m (6,<p) = (-l) m Y e * m (9, <p), and the facts Y em (9,0) G R and Y em (9,<p) = 
Y £m (9,0)e im v), 



Y lm {6,ip)Y lm {-,Q) (e- innpt 5 n , m -e inUpt 5 n ^ m ) 

n=-i 

I 

= -2i Yim{9, 0)Y em (|, o) sin(nfi P t - m^)5 n , m . (D.19) 



m=0 
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Then we are left with 

„ k A- 00 ^ 

e\x) = E E ^ - m ^ cos2 ^ Y *»( > (£> °) ( D - 2 °) 

x sin mfipr sm(mQpt — m(p)5 n ^ m . 

We now consider a regulated 0^. We truncate the sum over i at the cutoff L < £ crit . 
Since we are in the eccentricity case this will imply a cutoff in N because the kronecker 
delta forces n < I. 



VUJHpl * £=0 m=0 

x sin mQpr sin(mQpt — mip) . 



We are now interested in finding the maximum value of <p x > . This amounts to picking 
9 = 7r/2,<£> = and then considering points on the radial lightcone r = t with the 
specific values r = t = (2k + l)irQp/2, so that sin mQpr x sin mQpt = 1 for any even 
m. Call this set of parameters x max (note that x max is not a unique set of parameters, it 
is just any choice that satisfies these conditions). In this case (using the approximation 
for Y lm (it/2, 0)) 

d> { '(x max , £) ~ f= o = 7= L ■ (D.22) 

V 1 max ' ; 4.5 V6M P1 r2 7T 2 97rV6Mpir2 1 ' 

The factor of 1/4.5 is an approximation. Naively one would guess this factor is 1/4 
because we are taking roughly half of the m terms and half of the I terms. Numerically 
it seems that 4.5 for this factor is a better fit. 

D.3 Second Order Fluctuations 

For r > Op 1 we ignore the <ps piece compared to the 4>r piece, since we are interested 
in the derivatives of the field. We first compute explicitly the source for (p^ 2 \ 

J^ 2 \x) = A [ D7ro ( D 0(i)) 2 - nn (V M V,0«) 2 (D.23) 

-2V M V„7r o V At V> (1) D0 (1) + 2V M V iy 7roV i ,V A (1) V A V^ 
1 



A 3 



where the covariant derivatives are taken with respect to the flat 3d metric in spherical 
coordinates. As discussed earlier, since we work in the WKB regime we can do an 
expansion of derivatives in powers of 1/fip, d 2 f ~ ((nilp) 2 + ^^-f + •••)/■ Thus to 
leading order in Qp, 

/2, w-^(tr) 1/3 ^(-(^ (1, ) 2+ ^ ,11 )'>°- < DM > 
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and one needs to work to next to leading order in to see any relevant contributions. 
In this case the terms in the line of (D.23) that actually arose from C 3 are the dominant 
ones, 

J (2) = ~ [(<9 r 2 (1) ) 2 - 2 {d t d r ^) 2 + {dU {l) ) 2 } (D.25) 
= ^ [(W W ) 2 - (^ (1) ) 2 ] , (D.26) 

where we have used the fact that d r <p^ = d t (f>^ to first order in Qp 1 . For fixed 
£, £', m, vn! this looks like 

J ^^huw) Q " m m (D ' 27) 

TO IX 7T / 7T \ 

x cos 2 — Y £m (e : 0)Y im (-,0)Y ew (6,0)Y £ , m , (-,0j 

x ( cos mil pr cos m'VL pr cos(mf2 pt — rmp) cos(m'f2 pt — m'cp) 
— smmQpr sin m'VLpr sin {mVLpt — rmp) sm(m'Qpt — m'ip)) . 

Knowing the source, the equation of motion for the second order field fluctuation 

4>( 2 \x) is 

□0 (2) = (£) 2 (-^ 2 + 9 2 + i9 2 ) 0( 2 )(x) = / 2 )(x) . (D.28) 
In the WKB region r > Qj, 1 we expect the solution to take the from 

0( 2 ) rv, A(r, fi s ) sin(W + y? r ) sin(wt + ip t ) , (D.29) 

with A(r, fls) a slowly varying function. Considering the worst case scenario, we can 
neglect the angular dependence. In order to match with the source we take A(r) = ar 2 
leading to the ansatz 

2 L L e 

(2) 



* ' £=0 £'=0 m=0 m=0 

x (sin mflpr sin m'Qprsm(mQpt — rmp) sm(m'Qpt — rmp)) . 



To satisfy the equations of motion we must set 

2 / 2* Mr 



A 3 \V6M Pl r 2 



J npmm cos 2 ^Y £m (|, o) . (D.31) 



So now we have an expression for (p^ valid in the regime where £3 terms dominate. 
Now evaluating the second order field fluctuation at the set of point x max where the 
first order fluctuations take their maximal value, we get 

2 



£=0 m=0 



(D.32) 
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with the sum being approximated by 

L e r3 

E 2 ' 0) cos t * ' (D - 33) 

£=0 m=0 

so the second order field fluctuation at these set of points is 

As we see in section 3, (see eq. (3.35)), this usually implies a breaking of perturbation 
theory in many physical situations. 

E Validity of Perturbation Theory in the Cubic Galileon 

Finally we give a simple check on the validity of perturbation theory in the case of the 
cubic Galileon. Unlike in the quartic and quintic case, we find that the velocity pa- 
rameter v acts as a small parameter which ensures the convergence of the perturbative 
expansion. The properly normalized modes found in Ref. [35] are given by 

with 

f(M + l)/4fcr 
' [ -1/4 for f = K J 

We compute <pR in the WKB regime using the retarded propagator, as in appendix D. 
We find that for t > 

A/^A 1/2 ^ * ~ / r \ V4 cos (^nilpr + p) 

= -(w) E EE— («,)(- v ^ ; - 

\ / 1=0 m=—i n=0 x ' v 

(E.3) 

where P is an irrelevant phase factor. The time dependent coefficients a n e m (t) are 
given by 

• £ > modes 

arumit) = ^-n"' ( L ) 1/4 ^ x r^^{nn P {t-t'))e- imt 'f e {t'), (E.4) 
JWpi r * \ r */ Jo 2 p 

withM,, m = M 1 (M 2 /M tot ) (m)/2 +(-l) m M 2 (M 1 /M tot )^+ 1 )/ 2 , and a £ = Y, m (f,0)(^r + l)" 1 , 
and f e (t) = (1 - e 2 )(l + ecos(fipt)) -1 and by 

• £ = modes 

«noo(t) = j^~-n 7/A - x / ^ sin (nQ P (t - t')) f e {t') , (E.5) 
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where Ai is the reduced mass and where a' = F o(f , 0)2"( 7 / 4 )r(7/4)- 1 . The monopole 
and dipole must be treated separately from the other multipoles because at leading 
order in v the £ = and 1 part of 0# are zero by energy and momentum conservation. 
Thus the monopole and dipole both receive an extra power of v 2 suppression compared 
to the scaling v U£ from the other multipoles. 

Now we check the in the case of the cubic Galileon, the action for the perturbations 
4> contains, in addition to (2.19), the cubic term 

Sf = ~J d*x^- 3 Wfn<p, (E.6) 

which should be small in order for perturbation theory to be under control. To demon- 
strate this is consistent we consider the effect of these terms in the two distinct regions 
r > Qp 1 and r < Qp 1 . In the region r > flp 1 but r < r*(= r +3 ), a simple estimate of 
the ratio of the cubic to quadratic terms in the action gives, for £ > 1 

U Mr) ~ (u;f )^W /4 , (E.7) 



^(2) \r*J A 3 
which is small for all £ > 1 if we take r ~ ilj, 1 . For the monopole and dipole we get 



o(3) 



c(3) 

2 i D ^ 

~ v and 



1=0 



~ v A , (E.8) 



again evaluated at r ~ Vfp 1 in terms of the velocity v = VLpf. We see that the extra 
v 2 suppression for the monopole and dipole is crucial because it means that the ratio 
is also small when £ = 0, 1. Thus if we evaluate the flux on a sphere of radius r ~ Qp 1 
where linear theory is valid. 
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